Robust quantum information processing 
with techniques from liquid state NMR 

By Jonathan A. Jones 

Centre for Quantum Computation, Clarendon Laboratory, University of Oxford, 
Parks Road, Oxford 0X1 3PU, UK 

en' 
o . 

While Nuclear Magnetic Resonance (NMR) techniques are unlikely to lead to a large 
scale quantum computer they are well suited to investigating basic phenomena and 
developing new techniques. Indeed it is likely that many existing NMR techniques 
will find uses in quantum information processing. Here I describe how the composite 
rotation (composite pulse) method can be used to develop quantum logic gates 
■ which are robust against systematic errors. 
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1. Introduction 



Nuclear Magnetic Resonance (NMR) is arguably both the best technique and the 
worst technique currently known for implementing quantum information processing. 
The great strengths and weaknesses of NMR arise from the same fundamental cause: 
the low frequency of NMR transitions (typically around 500 MHz, corresponding 
Oh! to about 2 fj,eV). This makes NMR experiments easy to perform (the experimental 

timescale is conveniently slow), and also acts to minimise the effects of dccohcrence. 
However the extreme weakness of NMR signals means that it is not yet possible 
to study single nuclear spins: instead we must use macroscopic ensembles, which 
occupy hot thermal states. 

These strengths and weaknesses mean that, while it is extremely unlikely that 
liquid state NMR techniques will ever be used to construct a large scale general 
purpose quantum computer (Jones 2000), NMR provides an excellent technique for 
conducting preliminary studies (Cory et al. 1996, 1997; Gershenfeld & Chuang 1997; 
Jones & Mosca 1998; Chuang et al. 1998), and for developing techniques which will 
be used in large scale devices. Furthermore, the NMR community has developed 
a sophisticated library of techniques for manipulating nuclear spins (Ernst et al. 
1987; Freeman 1997; Claridge 1999), many of which can be directly transferred to 
manipulate qubits in other implementations. 

In this paper I discuss the method of composite rotations, also called composite 
pulses (Levitt 1986), which are widely used in NMR to combat systematic errors 
arising from inevitable experimental imperfections. While many composite pulses 
developed for use in NMR cannot be directly transferred to quantum computing 
some can be, and novel composite pulses have been developed specifically for use in 
quantum computing (Cummins & Jones 2000; Cummins et al. 2002). More recently 
the concept of composite rotations has been extended to two-qubit (controlled) 
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logic gates (Jones 2002). Although these robust quantum logic gates have been 
developed in the context of NMR, and are described here using NMR terminology, 
the basic ideas are entirely general and can be used in many other experimental 
implementations . 



2. Spins, qubits and the Bloch sphere 

The majority of conventional NMR studies are conducted on nuclei with spin I = 
1/2, and these also provide a natural method of implementing quantum information 
processing, as the two spin states, \a) and |/3), can be trivially mapped to the two 
basis states of a qubit, |0) and A general superposition state (Nielsen & Chuang 
2000) can be written as 

|V) = cos(6»/2)|0) + sin(e/2)e**|l) (2.1) 

(neglecting irrelevant global phases) and so can be thought of as a point on a 
unit sphere, traditionally called the Bloch sphere, with states |0) and |1) at the 
north and south poles, and the equally weighted superpositions lying around the 
equator. Any unitary operation on a single isolated qubit (any single qubit logic 
gate) corresponds to a rotation of this sphere. 

Single qubit rotations can be specified by their rotation axis and their rotation 
angle; the rotataion axis can itself be described by a single point on the sphere. 
Within the NMR literature spin states and unitary operations arc both described 
using the product operator notation (S0rensen et al. 1983; Ernst et al. 1987; Hore 
et al. 2000) . For a single spin 

1 ( ( cos 9 sin A e-^ 

2 V + \sin0e^ -cosfl 



MM 



1 • ( 2 - 2 ) 

= — (<t + sin 9 cos (f> a x + sin 9 sin <f>a y + cos 9 a z ) 

= \E + sin 9 cos (j)I x + sin 9 sin <f)I y + cos 9 I z 

where the product operators used in the final line are closely related to the corre- 
sponding Pauli matrices. NMR systems are usually in hot thermal ensembles, and 
so are not described by pure states but by highly mixed states (Jones 2001); for a 
single spin the thermal equilibrium state take the simple form 

P = 1 + eh (2.3) 

It is customary to neglect the 1 term (which cannot be observed by NMR tech- 
niques) and the e factor (which simply determines the signal intensity) and so the 
equilibrium state is described as I z . 



3. Pulses and pulse errors 

NMR experiments are composed of a series of radiofrcquency (RF) pulses, which 
cause rotations about axes in the xy-planc, and periods of free evolution, which 
for a single spin can be described as z rotations. Two particularly common pulses 
(Jones 2001) are the 90° pulse corresponding to excitation 

-i*/2Iy 

I z ► I x (3.1) 
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(closely related to the Hadamard gate) and the 180° pulse corresponding to inver- 
sion 

Iz ► -/* (3.2) 

(in effect, a NOT gate). In each case the rotation angle is determined by the applied 
RF power (which determines the rotation rate) and the length of time over which the 
RF is applied, while the phase of the rotation axis (in the xy plane) is determined 
by the phase of the RF field. 

Genuine experimental implementations are, of course, not quite as perfect as 
the description above implies. Clearly a rotation can go wrong in one of two ways: 
there can be errors in the rotation angle or in the rotation axis. The first type of 
error, usually called a pulse length error, typically arises when the RF field strength 
deviates from its assumed value, so that all rotations are systematically too long 
or too short by some constant fraction. Errors of the second kind, usually called 
off-resonance effects, occur when the RF field is not exactly in resonance with the 
transition, resulting in evolution around an effective field, tilted out of the xy plane 
towards the ±z axis. The unitary operation describing a real pulse therefore takes 
the form 

U = exp (-i x 2tt v x [(1 + g) (I x cos <f> + I y sin <j>) + fl z ] x t) (3.3) 

where v is the nominal rotation rate, t is the pulse length, (f> is the pulse phase, g is 
the fractional error in the RF power, and / is the off-resonance fraction, given by 
/ = 5/v, where 5 is the off-resonance frequency error. 

The ideal inversion pulse, equation (3.2), occurs when vt — tt, <p — 0, and 
/ — g — 0; in real pulses the last two conditions are relaxed. While both errors 
can (and do) occur simultaneously, one is typically dominant, and it is most useful 
to begin by considering these errors separately. Furthermore conventional NMR 
experiments begin in some well defined initial state, usually I z , and the quality of 
a pulse can be assessed by the overlap between the final state and its ideal form 
(for an inversion pulse, —I z ): it is not necessary (or even desirable) to consider the 
effects of the pulse on other initial states. 

In the presence of pulse-length errors, an inversion pulse performs the transfor- 
mation 

Iz > — cos{irg)I z + sm(ng)Iy (3-4) 

and the component of the final state along the —I z axis, cos(7rg), provides a con- 
venient quality measure, lying in the range ±1. This function is plotted in figure 1; 
clearly the sequence only performs well for very small values of g (near perfect 
pulses). The conventional NMR method for dealing with pulse length errors in in- 
version pulses is to replace the simple 180° pulse with the composite pulse sequence 
90° 180° 90°. This sequence has an inversion efficiency given by 

cos(7rg) + | sin 2 (7rg) (3.5) 

which is also plotted in figure 1. The composite sequence performs better than the 
naive sequence for essentially all values of g, but especially for moderate errors, 
in the range ±10%. The manner in which this improvement is achieved can be 
understood by examining the trajectory on the Bloch sphere in the presence of 
small errors (Freeman 1997; Claridge 1999). 
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Figure 1. Performance of naive (dashed line) and 90°180°90° composite (solid line) inver- 
sion pulses in the presence of pulse length errors. Performance is assessed by the component 
of the final state lying along the — I z axis. 

The situation in the presence of off-resonance effects is more complex, as both 
the rotation axis and rotation angle are affected. Composite pulse sequences are 
known which tackle these effects (Levitt 1986), but this point will not be explored 
further here. 

4. Pulses and logic gates 

Inversion, which takes I z to — I z is clearly closely related to the NOT gate, which 
takes |0) to |1), but the two processes are not simply equivalent (Jones 2001). The 
NOT gate corresponds to a 180° rotation around the I x axis, while inversion can 
be achieved by a 180° rotation around any axis in the xy plane. The difference 
is that an inversion sequence need only act correctly on the initial states |0) and 
|1) (corresponding to I z and —I z ), but a NOT gate must also act correctly on any 
superposition of these states. In NMR terms this means that the gate must also 
interchange I y and —I y , and must leave ±I X unchanged. It is, therefore, important 
to analyse the composite inversion sequence, 90°180°90°, to see how it performs 
with these initial states. 

In the absence of pulse length errors the composite pulse sequence does, in 
fact, perform correctly, but in the presence of errors the situation is not so good. 
While the composite pulse sequence performs better than a simple 180° pulse for 
±I Z states, it performs worse than the simple sequence for ±I X states; the effects 
of the two sequences on ±I y states are identical. This behaviour is exactly what 
one might expect: it seems intuitively reasonable that composite pulse sequences 
should redistribute errors over the Bloch sphere, rather than actually reduce them 
(Cummins 2001). If this were indeed the case, then composite pulses would have 
little to offer quantum information processing, but surprisingly some composite 
pulses arc known which perform well for all initial states. Such sequences, sometimes 
called Class A composite pulses (Levitt 1986), are of little use in conventional NMR, 
and so have received relatively little study. They are, however, ideally suited to 
implementing quantum logic gates. 

The first application of composite pulses to quantum information processing was 
by Cummins & Jones (2000), who used composite 90° pulses to reduce the influence 
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of off-resonance effects on an implementation of quantum counting. More recently 
Cummins et al. (2002) have described two families of composite pulse sequences 
which correct for pulse length errors. From here 1 shall concentrate on one of these, 
the BB1 sequence originally developed by Wimperis (1994). 



5. The BB1 composite pulse sequence 

The BB1 composite pulse sequence was developed with two principal aims: firstly 
to provide good compensation for pulse length errors, and secondly to provide a 
composite pulse sequence which could be used to replace any simple pulse at any 
position in a pulse sequence (Wimperis 1994). The second aim is essentially equiva- 
lent to seeking a Class A composite pulse, and BB1 does indeed have this property. 
The first aim is also well achieved by BB1, which provides a quite remarkable de- 
gree of compensation for pulse length errors: it is not only better than any other 
known Class A composite pulse, it can also provide better compensation that many 
conventional sequences tailored to specific operations (such as inversion). 

When assessing the quality of a Class A composite pulse it is necessary to deter- 
mine how well the unitary transformation actually implemented (V) approximates 
the desired unitary transformation (£7). A simple and convenient definition of this 
fidelity is given by 

|Tr(\/[/t)| 

T= Tr(UW) ^ 

(note that it is necessary to take the absolute value of the numerator as U and V 
could in principle differ by an irrelevant global phase shift). A simpler approach, 
appropriate to single qubit logic gates, is to note that any unitary operation on a 
single qubit is a rotation, and so can be represented by a quaternion 

q = {s,v} (5.2) 

where 

s = cos(0/2) (5.3) 
depends solely on the rotation angle, 9, and 

v = sin(0/2)a (5.4) 

depends on both the rotation angle, 9, and a unit vector along the rotation axis, a. 
The quaternion describing a composite pulse sequence is obtained by multiplying 
the quaternions for each pulse according to the rule 

qi * q2 = {si ■ S2 — Vi ■ v 2 , siv 2 + S2V1 + vi A v 2 } (5.5) 

while two quaternions can be compared using the quaternion fidelity (Levitt 1986) 

•^(qi, q2) = |qi • q2| = \si ■ s 2 + Vi • v 2 | (5.6) 

(it is necessary to take the absolute value, as the two quaternions {s,v} and 
{— s, — v} correspond to equivalent rotations, differing in their rotation angle by 
integer multiples of 2tt). For single qubit operations the two fidelity definitions 
(equations 5.1 and 5.6) are equivalent, and quaternions will be used from here on. 
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I shall take as my target operation a NOT gate, that is a 180° rotation; similar 
results can be obtained for any other desired rotation (Wimperis 1994; Cummins 
et al. 2002). Thus the quaternion representing the ideal operation is 

q = {0,(1, 0,0)} (5.7) 

while the quaternion representing the rotation which actually occurs (as a result of 
pulse length errors) is 

qi = {cos[(l + g)n/2], (sin[(l + g)n/2], 0, 0)} (5.8) 

giving rise to a quaternion fidelity of 

2 2 

Ti = cos(W2) » 1 - ^y- (5-9) 

(this expression neglects taking the absolute value, and so is only valid for values 
of g in the range ±100%). The conventional composite pulse sequence 90°180°90° 
which has the quaternion form 

q 2 = {sin 2 [ ff 7r/2], (cos[ff7r/2], -sin[ff7r]/2,0)} (5.10) 

gives exactly the same fidelity, JF 2 — cos(<?7r/2) = T\. This confirms that the con- 
ventional sequence does not actually correct for errors, but simply redistributes 
them (Cummins 2001). 

One BB1 version of a NOT gate takes the form 90^180^360^180^90^, where 
the phase angles fa and fa remain to be determined, and a phase angle of cor- 
responds with an x rotation. Note that this composite pulse sequence comprises a 
cluster of 180° and 360° pulses placed in the middle of a 180°, pulse. In the absence of 
errors the central cluster has no effect whatsoever, and the pulse sequence collapses 
to a simple 180°, pulse. In the presence of pulse length errors the central cluster will 
have some effect, and the intention is to choose values of fa and fa such that the 
effects of the central cluster compensate the errors in the outer pulses. Note that 
the sequence discussed here differs subtly from the original BB1 sequence described 
by Wimperis (1994) which had the cluster placed before the 180°, pulse; in fact it 
can be shown that the cluster may be placed at any point with respect to this pulse 
without affecting the fidelity (Cummins et al. 2002). 

The quaternion describing the BB1 composite pulse sequence is complicated. Its 
z component is zero, as expected for a time-symmetric pulse sequence (Cummins et 
al. 2002) , but the remaining components show a complex dependence on fa , fa and 
g. Progress is most easily made by expanding the quaternion as a Maclaurin series 
in g. The first order y component can be set to by choosing fa = 3 fa, leaving the 
approximate quaternion 

q 3 «{-Hl + 4cos^]2,(l,0,0)} ( 5 - U ) 

(neglecting terms 0(g 2 ) and higher). Finally the scalar part of q 3 can be made 
approximately equal to by choosing fa = ± arccos(— 1/4), and following previous 
practice the positive solution is taken. These choices result in a quaternion with a 
fidelity 

1 57T 6 6 

= 7777 (150cos(W2) - 25cos(3;77r/2) + 3 cos(5<?7r/2)) » 1 - -77777- (5.12) 
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Figure 2. Quaternion fidelity of naive (dashed line) and BB1 composite (solid line) 180^ 
pulses used to implement NOT gates in the presence of pulse length errors. The BB1 
sequence outperforms a simple pulse for all pulse length errors in the range ±100%, and 
is particularly good in the range ±25%. 



Table 1. Quaternion infidelities (1 — T) of naive and BB1 composite 180J pulses. 
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in which both the second and fourth order error terms have been cancelled. It is 
clear from this analysis that the BB1 composite pulse outperforms a simple 180° 
pulse for small values of g. In fact it does better for all values of g in the range 
±100%, as shown by the fidelity plots in figure 2; the most spectacular effects 
are seen within the range ±10%, as shown in table 1. BB1 pulses can implement 
extremely accurate gates in the presence of moderate errors (1-10%), while naive 
pulses require impossibly accurate control of the RF field strength to achieve the 
same quality. 

It is also instructive to examine the effect of the BB1 pulse on particular initial 
states. For initial states lying along any of the cardinal axes the BB1 sequence 
results in an error term of order g 6 , although the exact size of the term depends on 
the choice of axis. In comparison, for initial states along ±7 2 a simple 180° pulse 
results in an error of order g 2 , while the conventional composite pulse sequence, 
90°180°90°, gives an error of order g 4 ; thus the BB1 sequence acts as a better 
inversion sequence than the conventional composite pulse sequence designed to 
perform an inversion! For initial states along ±I y the conventional composite pulse 
sequence provides no compensation, and both it and and the simple pulse give errors 
of order g 2 . The only blemish on the BB1 sequence is seen when examining initial 
states along ±I X , for which the simple pulse performs perfectly (the conventional 
composite pulse gives an error of order g 2 ). This property of perfect behaviour along 
one single axis is a particular property of simple pulses, and cannot be achieved 
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with composite pulses. The very best behaviour for BB1 is observed for initial states 
along two particular axes in the xz plane, for which an error of order g 10 is seen. 

While the performance of BB1 is extremely impressive, it would obviously be 
desirable to find an even better sequence, with even better error tolerance. Although 
such sequences probably exist it is not clear how they can be found. Initial attempts 
in this direction (G. Llewellyn, unpublished results) have had no success, but have 
simply made clear how unusually good BB1 actually is. 

Very similar composite pulses can be obtained for other pulse angles (Wimperis 
1994; Cummins et al. 2002): a 9 a pulse is replaced by (#/2) o 18O^36O° 18O£(fl/2)o 
with <fi = arccos(— 9/An). There is, however, a subtle point concerning the accuracy 
with which such pulse sequences may be implemented. Typically all the pulses 
in such a sequence are implemented by applying the same RF field for different 
lengths of time, and the clock controlling the RF field has a finite time resolution. 
While it is not necessary to control the absolute lengths of each pulse to very great 
accuracy, it is essential that the relative lengths of each pulse are correct. This 
is easily achieved when 9 is 180° or some simple fraction of it, as all the pulses 
can then be implemented as multiples of some common element, but is much more 
difficult for arbitrary angles. 

6. Two qubit logic gates 

It is well known that any desired circuit can be constructed using single qubit logic 
gates in combination with any one non-trivial two qubit logic gate (Deutsch et al. 
1995). The two qubit gate most commonly discussed is the controllcd-NOT gate 
(Barenco et al. 1995) which applies a 180° rotation to its target qubit conditional 
on its control qubit being in the state |1). An essentially equivalent, and frequently 
more convenient, alternative is the controllcd-phase gate, which applies a 180° ro- 
tation to its target qubit conditional on the state of its control qubit. Note that 
in this case the logic gate acts symmetrically on the two qubits: the control/target 
distinction is convenient but artificial. This choice of two qubit gate is particularly 
convenient in implementations, such as NMR, built around Ising couplings, as the 
controllcd-phase gate and evolution under the Ising coupling are trivially related 
(Jones 2001, 2002). A controllcd-NOT gate can then be implemented by applying 
Hadamard gates to the target nucleus before and after the controlled-phase gate. 

Consider a system of two spin-1/2 nuclei, / and S. The Ising coupling gate is 
implemented by evolution under the J coupling Hamiltonian 

H IS =nJ2I z S z (6.1) 

for a time t = 4>/ttJ, where J is the coupling strength and <f> is the desired evolution 
angle. The desired controlled-phase gate requires cf) — n/2 and sot = 1/2 J; in 
NMR this is known as the antiphase condition. In order to implement accurate 
controllcd-phase gates it is clearly necessary to know J with corresponding accuracy. 
This is relatively simple in NMR studies of small molecules, but is much more 
difficult in larger systems. In particular, many experimental proposals contain an 
array of qubits coupled by Ising interactions (Ioffc et al. 1999; Cirac & Zoller 2000; 
Briegel & Raussendorf 2001; Raussendorf & Briegel 2001), with J couplings that are 
nominally identical but in fact differ from one another as a result of imperfections 
in the lattice. In systems of this kind it is desirable to be able to perform some 
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Figure 3. Pulse sequence for a BB1 based robust Ising gate to implement a con- 
trolled-phase gate. Boxes correspond to single qubit rotations with rotation angles of 
cj> = arccos( — 1/8) « 97.2° applied along the ±y axes as indicated; time periods corre- 
spond to evolution under the Ising coupling, ttJ2I z S z , for multiples of the time t — 1/4 J. 
The simple Ising gate corresponds to free evolution for a time 2t. 



accurately known Ising evolution over a range of values of J. Perhaps surprisingly, 
this is relatively easy to achieve using composite pulse techniques. 

The problem of performing accurate Ising evolutions is conceptually similar to 
that of correcting for pulse length errors in single qubit gates, and the solutions 
are closely related (Jones 2002). Ising coupling corresponds to rotation about the 
2I Z S Z axis, and errors in J correspond to errors in the rotation angle about this 
axis. These can be parameterised by the fractional error in the value of J: 

Jrcal - 1. (6.2) 



^nominal 



Errors of this kind can be overcome by rotating about a sequence of axes tilted 
from 2I Z S Z towards another axis, such as 2I Z S X . Defining 

6$ = exp[-t x 6 x {2I Z S Z cos0 + 2I Z S X sine/))} (6.3) 

allows the direct evolution sequence (7t/2)q to be replaced by the composite pulse 
sequence (7r/4)o(7r)^,(27r)3^(7r)^(7r/4)o with <f) = arccos(— 1/8). The tilted evolutions 
can be realised (Ernst et al. 1987) by sandwiching a 2I Z S Z rotation (free evolution 
under the Ising Hamiltonian) between (jj^ y pulses applied to spin S. After cancella- 
tion of extraneous pulses the final sequence takes the form shown in figure 3. Note 
that the labelling of the two spins as / and S is arbitrary, and the <j)^ y pulses can 
be applied to the other spin if this is more convenient. 

It is vital that any robust implementation of a quantum logic gate be built 
from components that are themselves robust. The robust Ising gate uses only two 
components: single qubit rotations around the ±y axes, for which robust versions 
are described above, and periods of evolution under the Ising coupling. As before it 
is not necessary to accurately control the absolute lengths of the five time periods, 
but they must have lengths in the integer ratios 1:4:8:4:1. 

The fidelity gain achieved for Ising coupling gates by this approach is identical 
to that achieved for single qubit rotations. (For two qubit gates it is necessary to 
use the propagator fidelity, equation 5.1, but as mentioned above this is equivalent 
to the quaternion fidelity in the single qubit case). As the controlled-phase gate 
corresponds to a 90° rotation, rather than the 180° rotations discussed in the case of 
single qubit gates, the fidelities are slightly different from those discussed previously. 
For a simple Ising gate 
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Figure 4. Propagator fidelity of simple (dashed line) and BB1 robust (solid line) Ising gates 
used to implement controlled-phase gates. The BB1 sequence outperforms a simple gate 
for errors in the range ±100%; over the range ±25% the robust gate is indistinguishable 
from perfection on this scale. 



while for a BB1 based robust Ising gate 



6„6 



637rV 
65536 



(6.5) 



The fidelities are plotted in figure 4. Clearly the robust Ising gate compensates 
well for small errors in J values, especially within the range ±25%. Over the range 
±10% the infidelity of the robust gate is always less that one part in 10 6 ; to achieve 
comparable fidelity with a simple gate it is necessary to control J to better than 
0.2%, more than 50 times more accurately then is needed for the robust gate. 



7. Conclusions 

Composite rotations show great promise as a method for combatting systematic 
errors in quantum logic gates. Without progress in this area attempts to build 
large scale quantum computing devices will founder on the need for impossibly 
precise experimental control. Methods have been derived for tackling both pulse 
length errors and off-resonance effects in single qubit gates (Cummins & Jones 2000; 
Cummins et al. 2002) and for tackling variations in the coupling strength in the 
Ising coupling two-qubit gate (Jones 2002). Together these provide a universal set 
of robust quantum logic gates within the Ising coupling model. Although developed 
within the context of NMR, Ising couplings play a major role in many proposed 
implementations of quantum information processing (Ioffe et al. 1999; Cirac & Zollcr 
2000; Bricgel & Raussendorf 2001; Raussendorf & Bricgel 2001), and these robust 
gates are likely to find their final applications elsewhere. 

I thank the Royal Society of London for financial support. I am grateful to S. Benjamin, 
H. K. Cummins, L. Hardy, G. Llewellyn, M. Mosca and A. M. Steane for helpful discus- 
sions. 
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